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We study the Casimir force between a perfectly conducting and an infinitely 
permeable plate with the radiation pressure approach. This method illus- 
trates how a repulsive force arises as a consequence of the redistribution of 
the vacuum-field modes corresponding to specific boundary conditions. We 
discuss also how the method of the zero-point radiation pressure follows from 
QED. 



* Revised summary of talk given to the ITAMP Topical Group on Casimir Forces ^ Cam- 
bridge, MA, March 16-27, 1998. 



'Understanding of signs is a sign of understanding. ' J. Sucher [1] 



At the end of his work on a force between two perfectly conducting parallel 
plates in vacuum as a consequence of the change in zero-point energy [2], 
Casimir proposed that this force can be interpreted as radiation pressure 
from the vacuum field. Later this interpretation was reintroduced by several 
other authors [3 - 5] and used [4, 5] for calculation of the Casimir force in 
such As was mentioned in [5], the radiation pressure approach can 

be systematically developed on the basis of QED. Actually this was already 
shown in 1969 by L. Brown and J. Maclay [6] (who also used this approach 
to compute the attractive Casimir force). We will discuss this question in 
some detail in Appendix. 

P. Milonni, R. Cook and M. Goggin [5] noticed a puzling character of the 
vacuum radiation pressure: In the case of two perfectly conducting plates the 
modes of the vacuum field in the space outside the plates form a continuum, 
corresponding to arbitrary wave vector k , whereas those inside are restricted 
to discrete values of kz (with the z-axis perpendicular to the plates). So 
there are "more"[| modes outside to push the plates together by the radiation 
pressure than there are modes between the plates to push them apart. This 
results in the attractive Casimir force^. However, they concluded that this 
argument is superficial, since it cannot explain a radially outward Casimir 
force, known in the case of a spherically conducting shell, in spite of the fact 
that there also should be "fewer" modes of the vacuum field inside the shell 

^Quotation marks since we compare infinite numbers! 
^This argument is also given in [7] and [8]. 
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than there are outside. 

In order to resolve the puzzle we should look more carefully to how bound- 
aries affect the zero-point field. Doing this one can notice that, loosely speak- 
ing, as a whole there are no fewer normal modes inside confined volumes than 
there are in in free space; they just are shifted to other frequencies. As it 
turns out this question was discussed in a rigorous way by G. Barton [9 - 10] 
and mentioned in [11]. Not knowing his results I came to the same conclusion 
from the following considerations. Consider, e. g., the space between two 
perfectly conducting plates separated by the distance I . As was mentioned 
above, the vacuum electromagnetic field is here a sum of modes (standing 
waves) permitted by the boundary conditions. Each such a mode is mathe- 
matically equivalent to a harmonic oscillator of the same frequency. When 
the distance between plates is very large we actually have the free space case. 
If now one adiabatically moves the walls toward each other, the "electromag- 
netic oscillators" will remain in the ground state; only their frequencies will 
be shifted to values corresponding to the changed boundary conditions [12]. 
Thus, the main effect of the boundaries is to redistribute normal modes : for 
some frequencies cu there are more modes than in free space, for others there 
are fewer. 

One can check this conclusion by comparison of the mode spectral densi- 
ties. In free space the mode spectral density in a volume V is 

while inside the cavity between two perfectly conducting plates [13], where 
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kz = mr/l, n = 1,2, 

pi(a;)=po(a;)^[l/2+E^(--H]- (2) 

^ m=l ^0 

Here c is vclosity of light, ujq = ctt/1, and 9 is the Heavisidc step function. 

In the case of a perfectly conducting plate parallel to an infinitely perme- 
able plate, with the same separation I, kz — {n+ l/2)n/l, n = 0, 1, 2, ... [14] 
and therefore the mode spectral density is 

p,{u;) = po(c^)- f:ei--m + 1/2). (3) 

One can see that for about a half of the frequencies and p2(<^) are 

greater than po{u!). 

It is just a result of the redistribution of normal modes that the pressure 
from the inside vacuum field Pout is different than the oppositely directed 
pressure from a free space vacuum field Pin ■ As noted above, for two con- 
ducting plates Pin > Pout- However, since the redistribution of modes de- 
pends on boundary conditions, the relation between two pressures in other 
cases can be opposite. 

As an example of such a case let us consider a perfectly conducting plate 
parallel to an infinitely permeable plate, mentioned above. If a plane wave 
has an angle of incidence 9 the radiation pressure exerted by such a wave on 
a plane, P — 2wcos^{9), where w is the energy density. So a vacuum field 
mode of frequency cu, which has an angle of incidence 9, makes a contribution 
to the pressure 

PM^2i^co.^m^^|(|f, (4) 
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where k = wc and V isa quantization volume. A factor 1/2 has been inserted 
because the zero-point energy of each mode is divided equally between waves 
propagating toward or away from each plate [5] . 

Therefore we find for the net pressure P — Pout — Pin, where Pout and Pin 
are vacuum radiation pressures directed outward and inward, correspond- 
ingly, the expression 

P = [he in^lY > / dK dky ' ^ / i 

' tJ^ "Jo \kl + kl + [{n + l/2)7:/lf)'^ 

poo poo poo ^2 

-(hc/TT^) / dk^ / dky / dkz r, (5) 

^ ' 'Jo "io 'Jo \kl + kl + kl)^ ^ ' 

which can be transformed using variables s = {l/7ry{kl+ky) and u = {l/7r)kz 
into the form 

°° roo rla 

n=o -^0 (s + (n-M/2)2)2 

-Tduu^r-^^]. (6) 

Jo Jo [s + u^) 2 

In order to regularize the divergent integrals in (6) a cutoff function 
fa{[s + u^Y^'^) must be introduced, with /„ — 1 when the parameter a tends 
to, say, zero: 

lim/„(p)^l, (7) 

a—>0 

herep = ([sH-m^]^/^). Further, one requires that fa{p) vanish rapidly enough 
for y — >■ oo , 



Jim/„(p)^0. (8) 



so that the function 



(s + m2) 



is finite and Ga{oo) = 0. We can thus rewrite (6) as 

POO 

P = Uma^o{hc7r/4l^)[Y^ Ga{n + 1/2) - / duGa{u)]. (10) 

Using the Euler-Maclaurin formula [15] one can find [12] that in the hmit 
a — > the difference in pressure is finite, independent of cutoff, and reduces 
to 

P = (7/8);ic7r/240Z^ (11) 

This coincides with the expression for the Casimir force for such a system 
obtained by Boyer [14] , who used the energy difference method and a special 
(exponential) form of cutoff function. 

Let us discuss shortly as to why the net vacuum radiation pressure has 
positive sign in the case under consideration. Let Ha{u) = Ga{n + 1/2) for 
n < u < n + 1 , where n is an integer. It follows from (10) that the net sign of 
the vacuum radiation pressure depends on whether the area under the step 
function Ha{u) is larger or smaller than the area under Ga{u). 

It is not difficult to see that if the curvature of Ga{u) is zero ior n < u < 
n + 1 , the area difference between Ha{u) and Ga{u) in this region is equal to 
zero. If the curvature is not zero the area difference will be bigger or smaller 
than these values depending on the sign of dP'Ga/dv? , i. e. on whether in 
the considered region Ga{u) is concave or convex. 

One can check that Ga{u) is primarily convex for any acceptable cutoff 
function /„. So the net area beneath Ha{u) is larger than that beneath 
Ga{u) and, therefore, in the case of one conducting and one permeable plate 
Pout — Pin > and we have repulsion. 
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So we showed in a very simple case how zero-point radiation pressure 
can lead to a repulsive rather than attractive Casimir force. It is the precise 
distribution of normal mode frequencies associated with specific boundary 
conditions, together with the fact that the function Fa{u) is primarily convex, 
that determines the sign of the force. 

Appendix 

Let us show how the method of the zero-point radiation pressure follows 
from QED and that it is equivalent to the method of the change in the energy 
of quantum fluctuations of the electromagnetic field. 

It follows from the operator Maxwell's equations that g, the operator for 
the momentum density of the electromagnetic field is 

g=i-(ExB-BxE), (12) 

where E = E(r,t) and B = B(r, t) are the Heisenberg operators of the 
electric and magnetic fields, g obeys the equation of continuity 

where T'-', i,j = 1,2,3 is the spatial part of the 4-dim. energy-momentum 
tensor operator T''^, //,A = 0,1,2,3 of the electromagnetic field, x° — ct, 
— x,x^ — y, x^ — z. T'*^ is a certain combination of electric and magnetic 
fields components and, therefore, is an operator function of x^, the form of 
which depends on the problem. We will consider below T'*^ renormalized in 
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such a way, that its expectation value in the ground state of the electromag- 
netic field, < T^^ >, is zero in free space.Q 

Eq.(13) describes the local conservation of the field momentum. Integrat- 
ing dg^/dt over a volume V and using the divergence theorem one can find 
the force exerted on this volume from the internal electromagnetic field. The 
expectation value for a force component, is 

dv<-^>=- dv< —— >= / < > rCda, (14) 
ot J ox^ Je 

where n is a unit outward normal to the surface which surrounds volume V 
and da is an element of the surface 

Let us show that determined by (14) coincides with the force deter- 
mined by the zero-point energy method: 

-*-|^' 

where 6e is a infinitesimal change in the zero-point energy of the electromag- 
netic field in the volume V under the infiuence of a virtual displacement (5a;* 
of its surface. 

Consider for simplicity the case of two parallel infinite plates separated 
by a distance I, with the 2;-axis perpendicular to them. Each plate is either 
a perfect conductor or is infinitely permeable. As follows from dimensional 
considerations, the energy per unite area of the electromagnetic field between 
the plates, e, is 

€(xl-\ (16) 



^The following is basically a slightly modified derivation given in [6]. 
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For a displacement of one of the plates along the 2;-axis 5x^ = Sz = 51. That 
is why the force per unit area acting on a plate in terms of the energy density 
w — e/l is 

= -| = 3^. (17) 
At the same time as follows from (14) 

Fz^ j <T^^> n^da, (18) 

where S is a surface of a plate, which has a unit area. It follows from the 
requirements T'^^ = , d^T^^ = and the symmetry of the problem that for 
the case under consideration [6] 

< T''^ Ci^g"^ - x^'x^) (19) 

where C is a constant, the metric g^'^ has the signature (—1,1,1,1), and 
x'^ = (0, 0, 0, 1) is a unit vector along the z axis. Thus, since g^^ — —1 and 



and 



w =< >= -^C, (20) 



r33^(l/3_i)^^_3^ (21) 



So, taking into account that < T^^^ is a constant we have from (19) 

F, =< >= —C. (22) 
Finally, from the comparison of this expression with (20) it follows that 

= 3w, (23) 
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which coincides with (17). 

In the case of two perfectly conducting or two infinitely permeable plates 
the renormalized energy density w is negative. So, as follows from (23), the 
force is attractive. In the case of a perfectly conducting plate parallel to an 
infinitely permeable plate w is positive and therefore the force is repulsive. 
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